Anomalous dielectric after-effect in ferroelectric KH2PO4 
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Dielectric permittivity, e = e' — if" of KH2PO4 pressed powders was measured between T = 1.4 
and 25 K in the presence of dc electric bias fields. Usually if the bias was switched at time t = 
0, e' and e" jumped to new values then decreased approximately as logt (6 s < t < 2000 s). This 
well-known effect, that is also found with single crystals, is attributed to switchable microdomains 
that are present accidentally in crystals but are systematically more numerous in powders. A very 
different after-effect was observed in a narrow T range around 7-8 K. t jumped to a lower value then 
increased with t according to a stretched exponential with a T-dependent time constant. This lay 
near the extrapolation of the Arrhenius law of a known, but unassigned weak-field dispersion that is 
a property of polydomain single crystals as well as powders. The weak-field dispersion is attributed 
to the elementary movement of a jog on a lateral step displacement of a domain wall, consisting of a 
single H-bond reversal. The anomalous after-effect results from the interaction between these point 
defects and the microdomain system. 
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I. INTRODUCTION 

It has been knownEI since the early days of ferroelectric 
KH2PO4 (KDP) that the material exhibits a retarded 
response to a disturbance as well as a prompt one. In 
particular, if its dielectric permittivity, e = e' — it" , is 
recorded using a low amplitude ac field in the additional 
presence of a dc bias field that is changed abruptly at time 
t = 0, then t' responds by promptly taking a new value, 
then by decreasing with logt as shown in Fig. |l|, curve 
(a). Zimmer, Engert, and Hegenbarthcl reported this for 
a single crystal with the fields parallel to the ferroelectric 
c-axis, at temperatures T = 4.2, 20.4, 77.8, and 300 K. 
The property is shared by many disordered ferroelectrics 
and will be referred to as "the normal after-effect" . 

Part of the interest of such studies lies in compari- 
son with dielectric glasses which also exhibit a "normal 
after-effect" at T < 1 KE It is also well known that ferro- 
electrics with diffuse transitions and relaxor ferroelectrics 
display low temperature thermal properties somewhat 
similar to structural glasses.Qcl A T 3 / 2 specific heat ter. 
at T < 5 K was also at one time reported for KDI 
though KDP has a normal, sharp ferroelectric transition. 
It was later clarified that large, pure KDP crystals can 
have almost Debye-like specificJaeatjS and lack a "glassy" 
thermal conductivity anomalyu 

The present article reports a study of the dielectric 
after-effect in KDP, mostly as pressed powders. These 
would be certain to have a strong T 3 / 2 specific heat term, 
and the after-effects were found to be stronger, and less 
variable from one sample to another than with single 
crystals. In the course of this study, a curious anomaly 
was found In a narrow temperature range around 7.5 K 
the effect of a bias switch was a prompt decrease of e' fol- 
lowed by an upward relaxation according to a stretched 
exponential law. Curves (b), (c) and (d) in Fig. I] show 
schematically what was observed at 9.0, 6.0 and 7.5 K 



respectively. The retarded response appeared to be the 
sum of two terms, the usual downward relaxation as logt 
(the normal after-effect) and an upward relaxation with 
a characteristic time constant that depended on T. Near 
7.5 K, the upward relaxation dominated throughout the 
range 6 s < t < 2000 s, at 9.0 K only at the shorter times, 
at 6.0 K the longer ones. This term will be referred to as 
"the anomalous after-effect" . Two electroded single crys- 
tal samples exhibited the anomalous after effect, as well 
as the normal one,0 but both effects were more precisely 
measurable using pressed powders. 

Among many previous dielectric studies of KDP, ones 
by Holste, Lawless and— SamaraEj and by Motegi, Ku- 
ramoto, and NakamuraOii 2 ] paid particular attention to 
temperatures below 25 K. Over a continuous background 
of absorption and dispersion rising regularly with T, 
Motegi et al. observed two specific dispersions, one sen- 
sitive to fields parallel to the ferroelectric axis {E || c), 
that required the presence of domain walls, the other 
sensitive to E \\ a and equally present in poly- or mono- 
domain crystals. The E || c dispersion (hereafter "KMN- 
C" ) obeyed an Arrhenius law with activation energy A |=j 
19 meV and pre-exponential frequency, /o = 45 GHz,E3 
i.e., it was centered at frequency / = fo exp(— A/fcT). 
It looks like the typical effect of a reorientable point de- 
fect species in a crystalline environment. A link between, 
the anomalous after effect and KMN-C was suggestedEI 
by the finding that the time constant of the stretched 
exponential obeyed the same Arrhenius law as KMN-C. 
As it will be shown in Sec. ^ below, the agreement is 
not as exact as in itially supposed, but it will also be re- 
ported (Sec. IV C) that the strengths of the two effects 
were closely correlated, when samples of different quali- 
ties were compared. 

The main aim of the present article is to report the 
anomalous after effect and to suggest an explanation for 
it. Since it appears to be related both to the normal af- 
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ter effect and to KMN-C, a report and discussion in some 
detail of each of these effects is also necessary. 

It is notorious that a dielectric study of any solid sam- 
ple that is not a properly electroded single crystal can 
lead to spurious results that bear little relation to bulk 
mate rial properties. It is argued below in Sec. Ill and 
IV C that serious errors are avoided by restricting atten- 
tion to T < 25 K, and interpreting the measurements 
cautiously. The general features of the data are also 
described in Sec. [II, by reference to one representative 
KDP sample. In Sec. IV A the normal after-effect is re- 



ported. This allows scaling parameters to be defined, 
that are useful in the subsequent (Sec. IV B and [V Oj ) de- 
scription of the anomalous after-effect. Results obtained 
with four single crystal samples are reported in Sec. 0, 
and the three effects, their interrelations and possible ori- 



gins are discussed in Sec. |VI 



II. SAMPLES AND MEASUREMENTS 

The pressed powder samples were derived from Aldrich 
99+% ACS reagent. In most cases the product was 
ground manually to ~ 3 fim particle size (as examined by 
swept-beam electron microscopy at 20 kV) then pressed 
(~ 50 MPa) into pills of diameter 10 mm and thickness 
normally 250-400 jum, density > 80% of crystal density. 
The pills were pressed (~ 100 kPa) between lead or in- 
dium electrodes. Four single crystal samples were also 
studied, and these are described in Sec. |V|, where their 
results are also reported. 

Capacitance and conductance were measured using a 
General Radio 1621 transformer bridge system, and ex- 
pressed as e' and e" based on the external dimensions of 
the pills and electrodes. The variations of e' and e" dur- 
ing a 360 s interval were obtained from a chart record- 
ing of the bridge off-balance signals. At longer times 
the bridge was rebalanced for each point. The tempera- 
ture was controlled using a carbon resistor sensor (Allen 
Bradley 390 Vt), that was calibrated periodically by sub- 
stituting a Pt and a Ge resistor for the sample capaci- 
tance. Thermometry errors were of two types, a short 
and a long term error. The short term error was caused 
by the carbon resistor varying with time following a tem- 
perature change. This was noticeable below 10 K and 
non-negligible below 6 K. It was allowed for by assuming 
that 
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where de' j dt was measured and (de' /dt)x was required. 
By waiting long enough before switching the bias field, 
the last term became a slow drift that could be subtracted 
confidently. The "long term" error had various contribut- 
ing causes. In routine work with different samples, the 
absolute T was known to ±0.5 K, but for the extensive 
work on Sample 2, to ±0.2 K. For the detailed study 



around 7.5 K (Fig. ||) the relative error was ±0.05 K 
(rapid succession of measurements without heating above 
122 K). For Fig. Ill absolute T was known to ±0.01 K. 



III. RESULTS, GENERALITIES 



The results reported here and in Sec. IV A were ob 



tained with a typical pressed powder sample (Sample 1) 
prepared from KDP as received (not recrystallized). 



A. Weak ac field, no bias 

Figure ^| shows the response of Sample 1 to alternat- 
ing fields of different strengths and a fixed frequency. At 
ambient T, sample impedance was limited by conduction, 
but this diminished rapidly on cooling, and became unde- 
tectably small already well above the ferroelectric transi- 
tion temperature, T c , and outside the range of Fig. ||. In 
Fig. ||, T c is marked by a peak of e' at the usual value of 
122 K. Well below T c there are e" peaks near 12 K and 
near 60 K. The 12 K peak is a genuine property of the 
material but the 60 K peak is not. A spurious peak is to 
be expected near 60 K for the following reason. For the 
sake of argument, suppose each grain had an anisotropic 
permittivity as measured with a properly electroded, un- 
damped singleperystal. Both e' c and e" c rise rapidly with T 
in this rangepEj e' c changing in order of magnitude from 
w 10 to > 10 4 . At e' c « 10, there would be appreciable, 
though nonuniform penetration of E \\ c field component 
into suitably oriented powder grains, but at e' c > 10 4 
such penetration would be negligible, the applied field 
being confined to the vacuum gaps, and to E || a within 
the grains (e' a remains moderate, rising to a peak value 

60 at T c , while e" remains low). The best compromise 
between the rising e"(T) and the falling E \\ c penetra- 
tion then locates a loss peak as in Fig. ||(a). It is not a 
Maxwell- Wagner effect and does not involve dc conduc- 
tion. 

On the other hand at T < 25 K, crystal e a and e c 
values are both moderate and depend weakly on T. At 
the moderate field values used, the dielectric response did 
not depend markedly on field strength (Fig. ||(a)), while 
dc conductance was totally insignificant. Insofar as the 
constituents of a composite dielectric material can be re- 
garded as continuous media, the e value of the composite 
is a mean value of the e of the constituent parts. The 
possible complex jneap values lie within rigorously de- 
termined boundsJlallS In these conditions, the 12 K loss 
peak in Fig. || can be assigned as a peak in e c or e a or 
both. In absolute value, it is distinctly stronger thap_ek 
ther the e c or e a dispersions reported by Motegi et alL3x3 
in single crystals, but in position it is nearer KMN-C (see 
Sec. V] below and Fig. |l2|) . The background value of e" in 
Fig. 2ns also an order of magnitude higher than the back- 
ground e" c in Motegi's crystals, so the relative strength of 
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the special absorption effect is quite similar. It is to be 
noted that the background e" value as well as KMN-C de- 
pended on the presence of domain walls )13 so the relative 
strength could be used for comparing the KMN disper- 
sion in samples of different qualities. It will be useful to 
characterize this relative strength for different samples 
by expressing the ratio of the peak e"(T) near 12 K to 
the minimum near 15 K, the ratio a/b in Fig. |^(b), as a 
percentage. For Sample 1, a/b = 126% and in absolute 
value, c « 0.042. 



B. Effects of bias switches 

Figure ^ shows typically how e' behaved following bias 
steps of different magnitudes at any T value not near 
7.5 K. "Step" here means a single abrupt change of ap- 
plied field. Normally after each step, and after the sub- 
sequent changes of e had been recorded the sample was 
heated to T c and recooled. As Fig. || shows, after a small 
bias step, e' promptly rose then gradually returned to- 
wards its original value, but after a larger step it grad- 
ually moved to new, lower values. A curve like (a) of 
Fig. |l| would be found after a bias step of intermediate 
magnitude. In each case, taking t = at the bias step, 
e'(t > 0) w e'(oo) + Ct~ p , where < p < 0.1 and C 
is a constant. To estimate e'(oo) would require a long, 
uncertain extrapolation, and over a limited t range the 
power law differs little from a logarithmic variation, so it 
is more useful to define s' — de' /dlni, and find s' , which 
was usually a slowly varying function of t. The variations 
of e" (not shown in Fig. |||) were similar, but scaled down 
by a factor of rj 4. Defining s" = de"/d\nt, \s"\ was 
like \s'\, usually a slowly diminishing function of t. In 
cases like curves (b) and (c) of Fig. [j], s' depended more 
strongly on t and even changed sign. In every case the 
final value of e', for t — ► oo, was less than or equal to the 
initial value e'(t < 0). The strength of the ac measuring 
field was unimportant provided it did not exceed 0.4A_B, 
otherwise \s'\ and \s"\ were underestimated. 

Figure || shows two relaxations of e'(t) that start from 
a same value, e'(t < 0) even though in one case the sam- 
ple had been cooled from T c in a bias field, the other 
in zero field. At T < 25 K, e'(t < 0) was generally re- 
produced to within 0.1% over a series of thermal cycles 
up to T c and back, even though the bias field was some- 
times zero, sometimes « 1 MV/m and the cooling rate 
also was variable. Similarly the value of e" prior to any 
bias step was reproduced to within 1%. This is quite dif- 
ferent from single crystal behavior. Three types of bias 
switches were studied. Either the sample was cooled in 
zero field, and field switched on at t = 0, or it was field 
cooled and at t = 0, the field switched off or reversed. 
Repeated measurements at 25 K with bias field switched 
from to 920 kV/m, separated by thermal cycles to T c 
yielded standard deviations ± 4% for s' and ±2% for s", 
excluding data from the first few cycles. If instead the 



sample was cooled in 920 kV/m and this bias switched 
off, \s'\ was 3±4% higher and |s"|, 2±2% higher. These 
are not significant differences, and the same applies if the 
sample was cooled in 460 kV/m and this switched to - 
460 kV/m. Only the magnitude \AE\ of the bias change 
was important. In this respect also the pressed powders 
behaved quite differently from single crystals. This point 
established, the three types of bias switches were used 
indifferently. Since also the sign of AE is irrelevant, it 
will be taken to be positive. 

A spurious after-effect could have been caused by space 
charge migration. If charge were to gradually migrate 
and accumulate, the effective penetration of the bias field 
into the grains would gradually diminish. This might 
cause e to return towards its original value, but not to 
acquire a new, lower value. It will be reported (Sec. IV C ) 
that samples of very different qualities always had lim- 
iting s' and s" values of the same order of magnitude, 
which argues against charge migration. 

Sequences of bias changes without change of T were 
not studied systematically, but schematically, if the bias 
field was switched periodically between values E\ and 
E 2 , starting at t = 0, without changing T, each switch 
initiated a new relaxation. In that case e'(t) rs e'(co) + 
C (t — ti)~ p , where ti is the time of the most recent bias 
switch. For small \E 2 — Ei\, e'(oo) w e'(t < 0), while for 
large \E 2 — Ei\, after several switches e'(oo) approached 
a constant lower value. 



IV. PERMITTIVITY CHANGES FOLLOWING A 
BIAS STEP 



A first systematic study of the after-effect at T w 25 K 
is reported because this is far from any special feature in 
Fig. ||(a), and the normal after-effect could be observed 
without the anomalous effect. 



A. Behavior near 25 K 

Figure ^ shows — s' and — s" vs. the magnitude AE 
of the bias step. Each pair of data points corresponds 
to a first bias step after a thermal cycle to T c and back 
to 25 K. s' and s" are the slopes of the best logarithmic 
fits to e'(t) curves as in Fig. || and corresponding e"(t) 
data for t between 6 and 360 s. It is already clear from 
Fig. that s' was not always proportional to AE. Fig- 
ure |4|shows that s' oc s" oc AE only at small values. At 
larger AE values s' and s" reached limits. The limiting 
value of —s' will be written s' , and another indepen- 
dent scaling parameter AE' will be defined by putting 
s'/s' = -AE/AE'q for small AE. and AE'^ simi- 
larly define the scale of the s"(AE) curve. The values 
of the four parameters for each of the three / values of 
Fig. [| are given in Table |. Their precise absolute val- 
ues are not significant on account of the incomplete and 
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nonuniform penetration of the applied bias fields into the 
grains of KDP. It is more useful to note that s' and s'q 
are decreasing functions of /, while AE' and AE' ' are 
increasing functions. Also, at given /, AE' 



1.5A££. 



TABLE I. Scaling parameters for the rate of change of per- 
mittivity of a typical pressed powder (Sample 1) at 25 K mea- 
sured between 6 and 360 s after a step, AE of bias field. / is 
the frequency of the low-amplitude measuring field. The rates 
of change are expressed as s' = tde /dt, s" = tde" /dt. Small 
steps, AE caused proportionate changes, s' / s' Q = —AE/AEq, 
s" /s'q = — AE/ AEq, but — s' and — s' ' represent limiting val- 
ues for s' and s". See also Fig. W. 



/ (kHz) 


AE' (kV/m) 
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AE ' (kV/m) 


II 

s 


0.12 


31 


0.0077 


20 


0.0018 


1.2 


37 


0.0055 


24 


0.0013 


12 


43 


0.0039 


28 


0.0009 



The curves in Fig. |J drawn to fit s' at 120 Hz and 
12 kHz correspond to an empirical formula s'/s' = 
-[l + {/AE' Q / AE) 2 ]- 1 / 2 with the s'q and AE' values given 
in Table [i]. The four other curves do not correspond to 
Table H, but were obtained from these two by supposing 
that for any given AE value, s' depends on the frequency 
of the measuring field according to a power law, and the 
response obeys the Kronig-Kramers relation. The ap- 
peal to Kronig-Kramers is based on the reasoning that 
whereas the response to AE on the time-scale of minutes 
is essentially nonlinear, e represents an approximately lin- 
ear response to a small field on the millisecond time-scale, 
and a set of e(/) data provides a "snapshot" for given AE 
and t. The same is true at t + St, and so of the derivative 
s' — is". It shows there is a link between the observa- 
tions that AE' Q and AEq increase with /, and that AE' Q ' 
is smaller than the corresponding AE' Q . 



B. T < 25 K 



The results reported here were obtained with another 
typical pressed powder (Sample 2) prepared from ma- 
terial that had been recrystallized in bidistilled water, 
the solution microfiltered. Referring to Fig. ||(b), a/b — 
130%, c w 0.028. 

Figure|| shows s'(AE) and s"(AE) for Sample 2 at a 
fixed /, fixed t interval and three temperatures. In this 
linear plot the details near the origin are not clearly seen, 
but the curves have two straight-line sections joined by a 
curved section. In a similar linear plot the same would be 
true of the 25 K data of Fig. ||, but at these lower T val- 
ues the straight-line section representing s'(AE) at high 
AE is not horizontal. It has a distinct positive slope at 
each T, steep at 7.5 K. In the case of s"(AE) the effect 
is much less pronounced. It is necessary to generalize 
the definit ions of the scaling parameters introduced in 
Sec. IV A and it will be defined that the straight lines 



intersect at (AE' ,— s' ) and (AEq,— s'q) respectively as 
shown in Fig. ^. 

The parameters so defined for Sample 2 at various T 
are shown in Fig. ^. s' and s' Q ' are increasing functions 
of T, but AE' and AE' ' peak near 5 K. Below 5 K, s'q 
and AEq, s' ' and AEq vary in the same proportions, be- 
cause the response to a small AE became temperature 
independent. 

In a linear plot like Fig. g, the slopes of the straight 
lines at high AE can be expressed in dimensionless units 
by using s' Q and AE' or s'q and AEq. In such units 
the initial slope is always -1 by definition. The values 
obtained for these positive slopes of s'(AE) were 0.01, 
0.02, 0.086,0.32, 0.04, and 0.017 respectively at T = 1.37, 
2.17, 4.9, 7.5, 9.9, and 12.4 K. For s"(AE) the slopes 
were 0.04, 0.09, 0.003, and 0.000 respectively at T = 4.9, 
7.5, 9.9, and 12.4 K. A closer scrutiny of the 25 K data 
showed a significant positive slope there also, for s' at 
AE > AEq. Statistical treatment of all data extending 
to AE > 20AE' o (five pressed powder samples) yielded a 
dimensionless slope 0.003±0.001. 

Figure || shows s'(T) and s"(T), at fixed / and fixed 
t interval, for Sample 2 and another similar sample at 
three different AE values. The AE values are such that 
AE > AE' (T) always, so that far away from 7.5 K, 
s' ~ — Sq and s" ~ — s'q. The position of the positive 
peak is independent of AE. In relative as well as in ab- 
solute value it is weaker in s" than in s' . 

Figure ^| shows s'(T) in the peak region at two frequen- 
cies and two time lapse intervals. The position depends 
distinctly on the time lapse but not at all on /. The mag- 
nitude is a very slowly diminishing function of /, which 
is consistent with the weakness of the s" peak (Kronig 
Kramers). 

Figure illustrates another method of studying the 
upward relaxation effect. Two data sets are shown, using 
different procedures. In one case all the e' measurements 
were made at 4.9 K, where the change Ae' caused by the 
bias step could be measured with little ambiguity be- 
cause s' was relatively very small. The sample was then 
annealed at 5.45 K for 2 mn and returned to 4.9 K for 
another e' measurement without further bias change, an- 
nealed at 6.0 K, at 6.6 K and so on. The x coordinate in 
Fig. |l^ is the highest anneal temperature prior to each 
measurement. The steepest slope near 7 K corresponds 
to the s' peaks in Figs. ^ and ^. Anneals up to 20 K con- 
tinued to have some effect, but then a plateau extended 
to 40 K, at which ~ 75% of the e' shift had been annealed 
out. The original e' value was almost entirely restored by 
an anneal to 65 K. For the other data set, e' was always 
measured at 9.9 K, where the bias step was effected, but 
after each 2 mn anneal at different T, and after return- 
ing to 9.9 K the sample was cycled to 122 K. Where they 
can be compared, ie T > 9.9 K the curves are similar 
apart from a scaling factor. The scaling factor suggests 
bias step AE was more effective at 9.9 K than at 4.9 K, 
doubtless because AE/AE' Q was larger. 

A similar procedure was used to compare the strength 
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of the anomalous after-effect of different samples after 
bias steps of different magnitudes. The bias was stepped 
at 6.0 K and the sample was annealed for 3 mn at 9.0 K 
before remeasuring e at 6.0 K. This anneal generally re- 
stored e' half way back to its original value. The Ae' 
recorded in these cases (and plotted in Fig. [ll]) was the 
change caused by the anneal alone, without reference to 
the original value before the bias step. 



C. Particle size and impurity effects 

The variations of the scaling parameters from one sam- 
ple to another are reported here, and also of the KMN 
dispersion strength and the strength of the upward relax- 
ation of e' corresponding to the anomalous after-effect. 
The principal results are also summarized in Table O. 



TABLE II. Summary of the effects of smaller particle size and of added impurities compared with "standard, pure" pressed 
powder samples. Different impurity species all had qualitatively similar effects. 



s 0i s o 

AE' , AEq 
KMN dispersion 

strength of anomalous after-effect 



Smaller particle size 

little changed (slightly increased) 

increased 

unchanged 

unchanged (at equivalent 
AE/AE' values) 



Added impurities 

unchanged 

increased 

weakened 

weakened 



Firstly it is necessary to consider variations amongst 
nominally pure samples prepared as in Sec. |ll|. This cate- 
gory includes Samples I and 2 that were said to be "typ- 
ical". It was found that s' Q , AE' Q , s'q, and AE' Q ' for such 
samples might vary by as much as a factor 2, but often 
much less. This would be due to accidental variations 
of density, homogeneity, and particle size distribution of 
the pressed powders. Therefore when comparing different 
categories of samples, any variations of these parameters 
exceeding a factor 2 are considered significant. Neither 
the scaling parameters nor the KMN strength depended 
significantly on whether the material was used as received 
or recrystallized, and whether in bidistilled water (solu- 
tion microfiltered) or in deionized water. There was no 
apparent difference between moderately pure and highly 
pure samples. Mean values at T = 25 K and / = 1.2 kHz 
were s'q « 0.005 and AEq w 35 kV/m. Referring to 
Fig. ||(b), the KMN dispersion strength was character- 
ized by I25%< a/b <I45% and 0.025< c <0.045. 

To investigate the effects of particle size, two samples 
were prepared from powders more thoroughly ground 
than usual, one of "as received" material, the other re- 
crystallized. Examination showed many particles of glob- 
ular shape and diameter ps I /xm. With the smaller aver- 
age particle size, these samples undoubtedly contained a 
higher proportion of severely damaged and nonferroelec- 
tric material. Both had slightly higher s' and s' ' values, 
but markedly higher AE' and AEq values than stan- 
dard samples. At T = 25 K and / = 1.2 kHz, s' « 0.007 
and AE ps 170 kV/m. On the other hand two sam- 
ples with larger than normal particles (loose powder con- 
tained angular shaped particles of dimensions « 10 ^m) 
yielded similar parameter values as standard samples. 
The larger particles would be likely to have broken up 
during pressing. Similarly, a normally ground sample 
pressed at 500 MPa had AE' and AE ' values typical of 
the more thoroughly ground samples pressed as usual at 



50 MPa. The KMN dispersion strength was found not to 
depend significantly on particle size. 

A series of samples was prepared in the usual way from 
material recrystallized from nonstoichiometric or impure 
solutions. Whatever the impurity species, the results 
were abnormally high AEq and AEq values, unchanged 
Sq or Sq and weakened KMN dispersions. In particu- 
lar, when the solution contained 0.2 H3PO4, 0.1 KHSO4, 
0.2 NH4H2PO4 or 2 RbH 2 P0 4 per 100 KDP, or 2 D 2 
per 98 H 2 (2% d), s' ranged from 0.0035 to 0.0077 as 
usual, but AE from 85 to 315 kV/m (always at T = 25 K 
and / = 1.2 kHz) while, referring to Fig. |(b), 101%< 
a/b <108% and 0.014< c <0.018. Increased impurity 
concentration (0.5 H3PO4, 0.3 KHSO4 or 11 RbH 2 P0 4 
per 100 KDP) caused further increase of AE' , still no 
change to s' but further weakening or disappearance of 
the KMN dispersion. Between 5 K and 25 K the tem- 
perature variations of the four parameters were roughly 
similar for pure or impure samples with coarse or fine 
grains (as Fig. [?]). 

The impure samples also had weaker anomalous after- 
effects. Figure [ll] demonstrates a correlation between 
strength of anomalous after-effect and KMN dispersion 
strength. Two groups of samples are featured. The first 
group comprises seven nominally pure samples. These 
were prepared as usual (Sec. |J) or were more thoroughly 
ground than usual so they had a variety of AE' val- 
ues, but all had normal KMN dispersions (a/b >125% 
and c > 0.025). The second group included the five 
impure samples mentioned above, and one other that 
fell into the same category (101%< a/b < 108% and 
0.014< c <0.018). As mentioned in Sec. |IVB|, the 



strength of the anomalous after-effect was characterized 
by Ae', the change in e' caused by anneal at 9.0 K, follow- 
ing a bias step at 6.0 K. To allow for variation of sample 
densities and homogeneities, Ae' was normalized with re- 
spect to s' (as measured at T « 25 K and / = 1.2 kHz). 
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If only the "normal KMN" samples are selected, Ae' /s' 
shows a strong linear correlation with AE, but only if 
the latter is also normalized with respect to AE . The 
line does not pass through the origin, but below it. For 
AE < 6AE' (AEq at 25 K, 1.2 kHz), or equivalently 
AE < 2AE (AEq at 7.5 K), the normal after-effect 
(downward relaxation) still dominated. The points in 
Fig. |ll| for the "weak KMN" samples fall near another 
line of lower slope, indicating a weaker anomalous after- 
effect. 

One other impure sample rdeserves a special comment. 
Following a known example ,113113 KDP was recrystallized 
from solution containing a large excess of base, 58 KOH 
for 100 KDP. The crystals were very hygroscopic and 
when removed from the dessicator and pressed into pills, 
fluid was expelled and filled the space between grains (pill 
density = 96% of crystal density). Several other samples 
were more or less conductive at room temperature, but 
in this respect, this one was an extreme case. Never- 
theless, it behaved at low T almost as a usual, slightly 
impure sample. At 25 K, e', e" , s' and s' ' were rather 
higher than usual. The ac field would have penetrated 
the grains more effectively, on account of the higher per- 
mittivity of the intergranular space. The KMN disper- 
sion was observed near 12 K (at 1.2 kHz) as usual, as also 
the anomalous after-effect centered around 7.5 K. This 
demonstrates conclusively than none of the reported low 
T effects is a purely surface effect, and that the KMN 
dispersion is unlikely to be a Maxatell- Wagner effect, as 
already argued by Kuramoto et aZJ13 

The 2% deuterated sample was classed amongst the 
impure samples because of its properties outlined above. 
It is logical to suppose that protonic impurity in KD2PO4 
(DKDP) would have a similar effect, so that if an ana- 
log to the KMN dispersion occurs in DKDP it would be 
necessary look for it in a d > 98% sample. Such an ef- 
fect was searched at 5 K< T <80 K but not found with 
a commercial (Aldrich) 98% d sample nor with another 
KDP sample recrystallized twice from 99.8% D2O in an 
atmosphere free of natural humidity. 



V. SINGLE CRYSTALS 

Four single-crystal samples were studied. SC 1 con- 
sisted of two slabs, cut normal to the c-axis and silver 
electroded by evaporation: total area 72 mm 2 , mean 
thickness 0.73 mm. The two pieces were connected in 
parallel. SC 2 consisted of two other c-cut slabs, gold 
electroded, also connected in parallel: 38 mm 2 x 0.35 mm. 
SC 3 idem but the two pieces were cut normal to an a- 
axis and silver painted: 75 mm 2 x 0.94 mm. SC 4 was a 
single c-cut plate of irregular shape 61 mm 2 x 0.30 mm, 
gold electroded. 

Both the normal (at 25 K) and the anomalous (near 
7.5 K) after-effects were observed with SC 1 and SC 2 
(E || c). All low T dielectric properties were sensitive to 



cooling speed through T c , and whether cooled in field or 
no field. They were less accurately reproducible from 
one thermal cycle to another than with pressed pow- 



ders. With SC 1, after fast zero- field cooling, 



and 



|s"| wereptypically several times smaller than with pressed 
powdersja but so also were e" and de'/dT. At 25 K, at- 
tempts to apply fields > 200 kV/m always caused an 
instability. At lesser AE values, \s'\ varied roughly as 
AE 0A , and no s' or AE' value could be estimated. Near 
7.5 K, the anomalous after-effect was characterized by a 
lower a constant in the stretched exponential law. Thus 
if e '(t) = e'(oo) - Cexp~(t/r) a , a took the value 0.45 
for pressed powder Sample 1 at 7.5 K, but 0.34 for SC 1 
at the same T. With SC 2, also fast zero-field cooled, all 
the low T dielectric properties were several times weaker 
than with SC 1. 

With SC 3 very little after-effect was found following a 
bias step, AE = 300 kV/m at 25 K. This puts an upper 
limit to \s'\ of 1.5 x 10~ 4 , so if any effect exists for E \\ a 
it is of a lower order of magnitude than for E \\ c. The 
inference is that the after-effects of the pressed powders 
are essentially caused by the E || c field component. 

SC 4 was used to check the Arrhenius law of the KMN- 
C absorption. For this purpose the calibration of the 
carbon resistor thermometer was not relied on, but a 
germanium resistor was placed inj-dose thermal contact 
with the sample. As KMN noted,li3 the background ab- 
sorption and dispersion that has to be subtracted de- 
pends on / as well as T. Some plausible choice has to 
be made, how this is done. KMN assumed at each T, 
a Cole-Cole law plus a linear e"(log/) background. For 
the present, the background to be subtracted at each T 
and / was assumed to be a linear interpolation of the 
data at T — 6.14 and 19.75 K (i.e. well below and above 
the KMN effect) and the same /. This also subtracted 
an instrumental and circuit error that becomes serious at 
/ > 10 kHz, and it resulted in a symmetric e"(log /). The 
e" peak position was found at each T. Two rfipld values 
were used, 3.5 and 7.0 kV/m (cf. 1.0 kV/m E3) and the 
data sets analysed separately. The 3.5 kV/m results are 
shown in Fig. 12 and yielded Arrhenius parameters fo — 
173 GHz, A = 20.09 meV. The 7.0 kV/m results would 
not be distinguishable in the plot and gave fo = 170 GHz, 
A = 20.03 meV. If the 6.14 K data alone had been sub- 
tracted as background, the figures would have been fo = 
605(640) GHz, A = 21.65(21.68) meV. Assuming the lin- 
ear interpolation is more appropriate, the presept result 



is not significantly different from the publishecEj result 
and it may be concluded that A = 19.5 ± 0.7 meV and 
40 GHz < f < 200 GHz. 

Also shown in Fig. [jj] are analogous results for pressed 
powder Sample 2, at 3.8 kV/m, using the same back- 
ground subtraction. It might be expected that e"(/) at 
6 K would have a peak near 100 Hz corresponding to the 
E || a dispersion)! 2 ! and so give a false baseline. Such 
an effect was searched between 5 and 7 K but not de- 
tected, which means it could have been no more than 
3% as strong as the absorption and dispersion near 12 K 
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in this sample. Best Arrhenius parameters for the latter 
were /o = 880 GHz, and A = 21.16 meV, and it is con- 
cluded from its position in Fig. [l^, and absence of E \\ a 
effect (at 5-7 K) that it corresponds almost entirely to 
the KMN-C effect. 

Figure [l^ also shows two points representing the time 
constant of the anomalous after-effect. The time con- 
stant t(T) was taken to equal t at the point of maximum 
positive slope of e'(logt). It was found at T = 7.44 and 
7.82 K, and (2ttt)^ 1 is plotted. The points fall distinctly 
below the Arrhenius law, whichever of the three data 
sets is extrapolated, but the law derived from the same 
pressed powder sample comes nearest. If an estimated 
normal after-effect contribution of the form Ct~ p had 
been subtracted from the data, the discrepancy would 
have been smaller. 



VI. DISCUSSION 

The after-effect in KDP at T < 25 K consists of two 
parts that are readily distinguished. Both, parts have 
been reported for single crystal samples, Bt!-as-. well as 
pressed powders and the KMN-C dispersionOllj is also 
a single-crystal property. The discussion of their origins 
must be broad enough to encompass both types of sam- 
ple. 



A. The normal after-effect 

A first salient feature of the normal after-effect is that 
it has been observed at T ranging from 1.4 to 300 K, 
so that its T range extends well above T c and well be- 
low the usual domain freezing temperature in crystals 
(~ 100 K). A second is that it has no characteristic re- 
laxation time, and a third is that its magnitude reaches 
a limit (s' , Sg) at very moderate values {AEq, AEq) of 
applied field step, with pressed powders at least. Below 
25 K, it conforms to a general law for "glassy" propex^ 
ties in that e is a simple, near-linear function of TlogtxB 
The low AE' and AEq values that decrease with T sug- 
gest a high dipole moment value that increases with T. 
All this suggests microdomains, that were present acci- 
dentally in the single crystals, but more systematically 
present in.-the powders. The notion of microdomains was 
developedcll to explain the diffuse nature of the phase 
transition in disordered ferroelectrics such as PLZT ce- 
ramics. TherSizes and shapes are widely distributed. R. 
Hohler et al£3 invoked this notion to explain their find- 
ing of a retarded dielectric response at low T in a PLZT 
ceramic. Although a different e was measured, related to 
the polarisation and not the etc polarisability, a T log t law 
was also found, between 20 and 80 K, in conditions cor- 
responding to s' = —s' Q (retarded response independent 
of AE). This possibly suggested independent thermally 
activated Debye processes with a uniform distribution of 



activation energies, which would imply a uniform spec- 
tral density of log r values (r is a relaxation time) that 
varies as T. Strongly interacting systems that relaxed ac- 
cording to a hierarchical sequenceE3 would also have been , 
possible, and more plausible. The authors suggestec£!rEij 
a discrimination in favour of the latter, based on results 
of applying a T step as well as a AE pulse, which in- 
dicated an abnormally low frequency prefactor. For the 
present, it is also more plausible to suppose strongly in- 
teracting systems in the pressed powder samples. This 
would explain why s' and Sg always took similar values 
at a given T. In the single crystals there may not have 
been enough active microdomains to reach the strong- 
interaction limit. Apparently AE' and AEq represent 
a threshold bias-step value for a prompt and widespread 
microdomain polarization rearrangement that limits the 
extent to which the system can be out of equilibrium. It 
is also most likely that microdomains are responsible for 
the "background" dispersion and absorption at low T, 
both in single crystals (E || c) and powders. They also 
explain naturally why AE' Q and AEq decreased with T, 
but increased with /. Either increased T or decreased / 
would bring into play larger microdomains. 

The microdomains that are active at low T are proba- 
bly ones that have less than full orthorhombic distortion, 
so they are easily switched. In single crystals they are 
likely to be associated with crystal defects, in pressed 
powders likely to be much influenced by the damaged 
grain surface layers. In the powders they must be sub- 
ject to a very broad distribution of stresses, as witnessed 
by the equivalence of field-cooled and zero-field-cooled 
properties. On the other hand, the active microdomains 
at ambient T are perhaps to be identified with the or- 
thorhombic inclusions that have been reported in the 
tetragonal phasc.Eil 

A "normal" after-effect of a bias step is also character- 
istic of disordered dielectrics generally at very low T. The 
immediate rise of e' and e" followed by decay as logi was 
observed with hydroxyl-doped KC1 at T <■ 1 KfB and 
several structural glasses, also at T < 1 KB The results 
could be displayed in plots like the present Figs. ||and||, 
though the samples were not thermally cycled between 
measurements as in the present study. The behaviour 
of the structural glasses below 1 K was explainedc3 by 
reference to a random Ising model of dipoles with long- 
range interactions .Ell One difference from the present re- 
sults was that the AE' a parameter was an increasing func- 
tion of T. In accordance with the model, it was found 
AEq w kT/p, where p is the (fixed) relevant dipole mo- 
ment. Another difference was that s' Q decreased sharply 
with T. The other main point in common is that e' and 
e" relaxed downwards, and this can be understood as the 
system of interacting dipoles gradually self-trapping into 
deeper potential wells, from which it can respond less 
actively to weak applied fields. 
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B. The anomalous after-effect 



C. Possible origin of KMN-C and anomalous 
after-effects 



The most obvious distinguishing features are that e' re- 
laxed upwards, not down and that although some effect 
was detectable at all T between 1.4 and 25 K, there was 
a pronounced peak near 7 or 8 K, depending on the time 
lapse. Also, the anomalous effect was only observed when 
AE > AEq, and then its strength as represented for ex- 
ample in Fig. [ll] increased as (AE / AEq — a), where awl. 
The AE' parameter needed to describe the normal after- 
effect appears as a threshold for the anomalous effect. 
More precisely there was a tendency for the strength to 
level off and possibly .saturate at some high AE value not 
reached in this work.la This suggests that the unit dipole 
moment involved here is much smaller than the moments 
of the microdomains. Since the anomalous after effect is 
also linked to KMN-C, it may be supposed a same species 
of dipole ("the KMN dipoles") is responsible. 

A phenomenological model can then be formulated as 
follows. Two subsystems (the microdomains and the 
KMN dipoles) coexist and interact, each having its own 
relaxation dynamics. The KMN dipoles only interact 
weakly with one another, but certainly interact with the 
microdomain system. After a bias step that is big enough 
to unsettle the microdomain system, this settles into 
a new configuration that is metastable, with the KMN 
dipoles in their present states. These dipoles then relax 
with their characteristic time constant, r, perpetually 
changing the local strain fields and electric fields. This 
affects the microdomain system qualitatively as would 
a random series of applied field changes. Referring to 
the effe ct of repeated bias changes that was mentioned 
in Sec. IIIB, and summing over the sample volume, it 



might be expected that after a single applied field step, 
AE > AE' , 

e'(t) = e'(oo) + Ce~ t/T t- p + — f* e - {t - u)/T (t - U)~ p dU 

T JO 

For simplicity, unstrctched exponentials have been writ- 
ten here. The term in C represents regions of the sample 
where e'(t) is still relaxing downwards due to the applied 
field step. Its volume diminishes exponentially. The C 
term represents the sum of micro-regions that have been 
affected by a subsequent KMN dipole flip, that occurred 
at ij. The expression allows an upward relaxation of e'(t) 
only if C > C. This is possible because if AE > AE' , C 
takes its limiting value related to Sq, while C can plau- 
sibly exceed this limit. The field changes caused by the 
dipole flips may be very strong, but they are localized and 
random, so they do not cause a prompt and widespread 
microdomain rearrangement like a strong applied field 
change. The microdomain system can therefore be driven 
further away from equilibrium. 
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KMN suggestecO that a peculiar mode of motion re- 
lated to the domain wall structure may be responsible. 
If some of the atoms located at the domain boundary 
were in shallow double potential wells, movement of these 
atoms across the barriers would be a possible origin. This 
brings to mind the hydrogen atoms and their O. . . H-0 
bonds. At the time of KMN, the potential barrier.-tp H- 
bond reversal was thought to be well over 100 meV,c3 but 
newer workE3~EHI based on neutron incoherent scattering 
data has yielded the much lower estimate of 37.1 meV,El 
and again a precise potential function with a barrier 
height near 30 meV (for DKDP, 135 meV) has been cal- 
culated for the system comprising a hydrexgen atom and 
coupled lattice mode, based on 30 K data.ES The princi- 
pal aim was to account for T c and its deuteration shift 
by considering lattice dynamics at T > T c , but the result 
might also apply to a special case where a single H bond 
reversal could occur at low T without the energy cost of 
creating a HPO4 - ion and an adjacent H3PO4 group. 

It is necessary to examine possible domain-wall struc- 
tures. Wall width increases with T towards T c ,Ej but at 
T C T c , walls are plausibly approximated by the va: 
ishingly thin models proposed by Barkla and Finlayson 
There were two such structures, the "polarized" and the 
"neutral" domain wall. Both respect the Slater rule (no 
HPO4 and no H3PO4 groups). It is not known for certain 
which is the more accurate approximation to real walls 
at T <C T c but Bjorkstam and Oettel calculated that tha 
polarized wall would be the more stable.La BornarelE3 
showed that either of these walls, normally planar and 
perpendicular to an a-axis, can have any number of lat- 
eral step displacements. The minimum displacement is a 
half lattice parameter. The step displacements, or quasi- 
dislocations, like the planar walls, respect the Slater rule, 
but only if they run straight across the entire crystal, par- 
allel to the c-axis. They are associated with an intense 
local strain field. Minimum domain wall movement there- 
fore involves glide of a quasi-dislocation along its entire 
length, or the presence of HPO4 and H3PO4 groups. Oth- 
erwise, if the walls have finite width, a minimum move- 
ment without HPO4 or l^PCL-groups consists of six si- 
multaneous H-bond reversals.EJ The simplest case with 
HPO4 or H3PO4 is a unit jog, where a quasidislocation 
shifts by one lattice parameter. This requires a single 
HPO4 or H3PO4 group and its illustration^ is repro- 
duced in Fig. [l3] for the case of a polarized domain wall. 
The case of a neutral wall is almost entirely equivalent. 
By a sequence of single H-bond reversals, the vacancy (or 
the excess H) can move from one PO4 to the next within 
a c-stack of unit cells, and so, in principle right across 
the crystal, together with the associated jog. In practice 
it may get pinned somewhere along the line, by a crys- 
tal dislocation or an impurity. A few such pinned jogs 
may persist as the material is cooled down to tempera- 
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tures where their creation by thermal activation would be 
virtually impossible. Also, progression of a jog requires 
successive reversals of differently oriented H-bonds, first 
one that lies near the plane of the domain wall, then 
one nearly perpendicular to it and so on. It is plausible 
that sometimes, because of the local fields and strains 
one of these, but not the other remains possible at low 
T, because in just one case the two states are near en- 
ergetically equivalent. This would constitute the double 
potential well system— 

Yamada and IkedaEJ considered the "cluster tunneling 
mode" or protonic polarontJ model, but found that their 
purpose of predicting hydrogen dynamics at T > T c was 
better served by a model of incoherent tunneling between 
self-trapped states. At low T this would become coherent 
phonon- assisted tunneling, with an expected transition 
probability oc T 7 . On the other hand, the protonic po- 
laron would have an extremely small tunneling splitting, 
so that thermal activation down to 12 K, or even 7 K, 
would be more plausible. Moreover, the tunneling mode 
would be overdamped at higher T, and in that condition 
the predictions of the two models would be experimen- 
tally indistinguishable. 

In anoiher development since KMN, evidence has been 
reportedE3 of an orthorhombic-monoclinic phase transi- 
tion in KDP near 60 K. This would mean the H-bonds 
within a ferroelectric domain are not all equivalent, as 
previously supposed, but not enough is known about the 
new phase to draw any other conclusions. 

Meanwhile, a different origin for the KMN-C effect can- 
not be totally rejected. This would attribute it to a defect 
species that is intrinsic to solution grown KDP (for exam- 
ple a growth dislocation or an included water molecule). 
At low T this defect would only be dielectrically activated 
by the presence of a domain wall. 



VII. CONCLUSIONS 



by point defects in crystals, and are unknown in struc- 
tural glasses at low temperatures. Before the present 
work there was no apparent link between KMN-C and 
the normal after-effect. The KMN-C dispersion possibly 
owes its origin to rare, isolated HPO4 and H3PO4 groups 
associated with jogs on lateral steps (quasidislocations) 
of domain walls. The activation energy of 19 meV, which 
is also shared by the anomalous after-effect, would then 
be related to the energy barrier for the reversal of a sin- 
gle H-bond in this particular environment. It is not clear 
how closely this should be assimilated to the hypothetical 
barrier for the reversal of a single H-bond in tetrag 
KDP, that has been the object of recent estimates! 
A correlation might be expected, if different ferroelectric 
compounds of the KDP type could be compared. The 
most interesting compariso n wou ld be with DKDP, but 
so far, as mentioned in Sec. IV C no dispersion analogous 
to the KMN effect has been observed with DKDP. For 
the other isostructural compounds the situation can be 
summarized as follows. As pressed powders, RbH 2 P04 
and the arsenates all behaved analogously to KDP, but 
the arsenates had first to be crystallized with excess base 
(solution pH>6). Each compound exhibited a low-field 
dispersion that obeyed an Arrhenius law, and a corre- 
sponding anomalous after-effect. With RbH2P04 the ac- 
tivation energy was close to the value for KDP, perhaps 
1 meV higher. With KH2ASO4 it was near 30 meV, with 
CSH2ASO4, 44 meV and with RKH^AsO^ between these 
two. However, no independent estimates of the barrier 
heights are available. 

The anomalous after-effect results from an interaction 
between two coexisting subsystems. One subsystem re- 
laxes with a definite relaxation time, the other with a 
very broad distribution of relaxation times. The result 
of the interaction is an apparent tendency of the system 
to evolve temporarily away from its stable equilibrium. 
The lines of an explanation sketched in Sec. VI B need to 
be developed into a model. 



A new dielectric property of ferroelectric KDP has 
been reported, "the anomalous after-effect" . It-has been 
shown to be related to the "|BO¥*nal" after-effectB and also 
to the KMN-C dispersion.llilll3 The normal after-effect 
is a well-known property of c-cut KDP crystals, and of 
many other ferroelectric materials. Like the T 3 / 2 specific 
heat term in microcrystalline KDP, it is reminiscent of 
such an effect in structural and dipole glasses, but differs 
in certain important details. It is attributable to mi- 
crodomains with a wide distribution of sizes, shapes and 
stresses. Study of the normal after-effect with pressed 
powders has allowed a set of scaling parameters to be de- 
termined for a series of samples of different qualities. As 
shown in Fig. [ll], these scaling parameters also apply to 
the anomalous after-effect, which therefore also involves 
microdomains. 

The KMN-C dispersion is another known property of 
cut KDP crystals. Such dispersions are generally caused 
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FIG. 1. Schematic illustration of dielectric after-effects. 
The field applied to the sample consists of an ac measuring 
field of fixed low amplitude plus a dc bias field. After the 
sample has been stabilised at a given temperature for some 
time, the bias field is switched. The real (in-phase) part of the 
permittivity, e', changes abruptly then usually relaxes down- 
wards (a). It may also relax upwards at short times (b), at 
long times (c), or at all times accessible to experiment (d). 
A slope s' — de'/dlnt is defined as the best fit to the data 
within a specified range t\ < t < t2- This is particularly 
useful for cases like (a) where s' takes a negative value that 
only depends weakly on t\ and ti- The imaginary (quadra- 
ture) permittivity behaves similarly but on a reduced scale, 
and s" = de" /dint. 
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FIG. 2. (a) Real and imaginary dielectric permittivity of 
a typical pressed powder sample (Sample 1) at 1.2 kHz and 
three field levels. Here and in Figs, ^j, ^| (?], and ^| e or related 
data are shown as filled symbols, e" data as corresponding 
outline symbols. The present data were not sensitive to tem- 
perature history on a scale that would be visible on the plot, 
provided the sample was held at each T value for 10 minutes 
or more, e is shown only for 4.3 kV/m as the other data 
would be indistinguishable on the plot, (b) Low T detail (e" 
at 1.4 kV/m only). The strength of the specific absorption, 
c relative to the background might be expressed by the ratio 
c/d, which requires an interpolation of the background curve. 
For comparison of different samples the ratio a/b is used. 
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FIG. 3. Time dependence of real permittivity of a typical 
pressed powder (Sample 1) at 25 K measured using 4.3 kV/m 
at 1.2 kHz. Two sets of data are shown. Points A: the sample 
was first cooled from above 122 K to 25 K in zero bias field 
and held there for 1 hr, then at t = a bias field of 11.4 kV/m 
was switched on. Points T: the same, except that the sample 
was cooled in a bias field of 400 kV/m and this field was re- 
versed at t = 0. e'(t) depended solely on the magnitude (11.4 
or 800 kV/m) of the bias step, and not on the value (zero or 
nonzero) of the cooling field. 
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FIG. 4. Sample 1 at 25 K in weak ac fields: derivatives of 
real and imaginary permittivity with respect to Int where t 
is the time since a bias step of magnitude AE. The sample 
had been cooled from above 122 K and held at 25 K for at 
least 1200 s when the bias was stepped. Data points are best 
fits over the interval 6 s< t < 360 s. T, •, A (for e data), 
V, o, A (for e" data) respectively / = 120 Hz, 1.2 kHz, and 
12 kHz. All de /dint and de" /dint values were negative, as 
in curve ( a) of Fig. hi. The fitting curves are defined in the 



bias step, AE (MV/m) 

FIG. 5. Data for a typical pressed powder (Sample 2) in 
weak 1.2 kHz field: derivatives of real and imaginary permit- 
tivity with respect to In t where t is the time since a bias step 
of magnitude AE. Data ponts are best fits over the interval 
10 s< t <100 s. T, », A (for e'), V, o, A (for e") respectively, 
T = 4.9, 7.5, and 12.4 K. Note that unlike Fig. §both scales 
are linear and that the scale of de" /dint is expanded and 
shifted. 
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FIG. 6. Schematic representation of data as in Fig. |H] to 
show the general definitions of the scaling parameters. In case 
of a "normal after-effect" with no "anomalous after-effect", 
the rectilinear parts at high AE would be horizontal. 
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FIG. 7. Scaling parameters as defined in Fig. ^, for Sam- 
ple 2, using a weak 1.2 kHz measuring field. The parameters 
are primarily related to the normal after-effect, but the val- 
ues obtained for s' near 7.5 K may have been influenced by 
the anomalous effect. Note that for disordered dielectrics-.at 
T < 1 K, AEq tx T, while s' sharply decreased with T.0O 
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FIG. 9. Similar plot to Fig. H showing details of the anoma- 
lous positive peak. All data here are for Sample 2 in 3.8 kV/m 
fields of two different frequencies, / following bias steps of 
magnitude AE = 2.25 MV/m. / = 39 Hz (□, 0), 10 kHz (■, 
Data points are best fits over the interval 10 s< t < 35 s 
(■, □), or 100 s< t < 330 s (♦, 0). de'/dlnt reaches a peak 
at a temperature that depends on t, not on /. 
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FIG. 8. Derivatives of real and imaginary permittivity (in 
weak 1.2 kHz fields) with respect to lnt where t is the time 
since a bias step of fixed magnitude, AE, applied at different 
T. Data points are best fits over the interval 10 s< t < 100 s. 
• , A, T (for e'), o, A, V (for e") respectively Sample 2 with 
AE = 2.25 MV/m, another similar sample with AE = 1.32, 
0.66 MV/m. Note that the scale of de" /dint is expanded and 
shifted, de'/dlnt has a negative term, that varies roughly as 
T ("the normal after-effect") and a positive term that peaks 
at T « 7.5 K ("the anomalous after-effect"). 
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FIG. 10. Anneal-out of the change in real permittivity of 
Sample 2 (in weak 1.2 kHz field) caused by a bias step of 
2.25 MV/m at 4.9 K (♦) or 9.9 K (■). Di fferen t procedures 



were used as explained in the text (Sec. IV B) but e was 
always measured at 4.9 K or 9.9 K, never at the anneal tem- 
perature and Ae' is the change in e' resulting from the bias 
step and any subsequent annealing. The anneal temperature 
is either the temperature of a single 2 mn anneal (■) or the 
highest, and latest of a series of 2 mn anneals (♦). 
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FIG. 11. Strength of anomalous after-effect vs. bias step 
magnitude for thirteen pressed powder samples. The bias 
was stepped at 6.0 K, the sample was annealed for 180 s at 
9.0 K then returned to 6.0 K. Ae' is the change in real per- 
mittivity (weak 1.2 kHz field, at 6.0 K) caused by the anneal. 
Where the normal after-effect was dominant it is negative, 
where the anomalous effect dominated, positive. Ae' is plot- 
ted normalized with respect to the scaling parameter s' (at 
25 K, 1.2 kHz). The samples form two groups: pure sam- 
ples with KMN dispersions of normal strength (•), impure 
samples with weak KMN dispersions (+). In the main plot, 
the bias step magnitude is also normalized with respect to 
AE' (25 K, 1.2 kHz). If it is not (inset, points • only) no 
clear pattern emerges. This is a key result. The anomalous 
after-effect is seen to correlate with the KMN dispersion, but 
also to be linked to the normal after-effect since the same 
scaling parameter, AE' is involved. 




PO4 tetrahedra and their hydrogen atoms 



TsL' 



ConventionalB'EHEl 



polarisation 
into page 



FIG. 13. (a) ConventionaltJ'OTJ representation of 
KH2PO4 showing a polarized domain wall with a lateral step 
displacement. Black dots represent H atoms. The dashed line 
encloses a stack of (tetragonal) unit cells that lies astride the 
domain boundary. Within this stack, the PO4 tetrahedra lie 
on a helix, whose pitch is the lattice parameter, (b) the same 
stack viewed at a higher level. To pass from one arrangement 
to the other, any one of the tetrahedra within the stack must 
possess only one closely bound H (the ion is HPO^ - ). The 
fault can travel up or down the helix bji-a series of single 
H-bond reversals as explained by BornarelEZI An opposite jog 
would require one H3PO4, and would be similarly mobile. 
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FIG. 12. Arrhenius plot showing data from Single Crystal 
Sam»l« SC 4, Pressed Powder Sample 2 and from Kuramoto 
et alL3 o, frequency of absorption peak of SC 4 at 3.5 kV/m, 
with background subtracted as stated in text; +, idem of 
Pressed Powder Sample 2 at 3.8 kV/m; •, (2ttt)^ 1 of Pressed 
Powder Sample 2, where r is the value of t at the inflexion 
point (maximum positive slope) on the plot of e' vs. \ogt 
after a bias step of 2.25 MV/hi. The line labeled C is the 
Arrhenius fit given by KMN,E3 and the curve labeled A an 
approximate fit to their E || a data. 
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